Mutual information and Bose-Einstein 
condensation 



C. N. Gagatsos*^ A. I. Karanikas''"^, and G. I. Kordas^^'^ 

^Quantum Information and Communication, Ecole Polytechnique de 
Bruxelles, Universite Libre de Bruxelles, 1050 Brussels, Belgium 
^ University of Athens, Physics Department, Nuclear & Particle 
Physics, Section Panepistimiopolis, Ilissia GR-15771, Athens, Greece 
^Institut fiir theoretische Physik and Center for Quantum, Dynamics, 
Universitdt Heidelberg, D-69120 Heidelberg, Germany 



Abstract 

In the present work we are studying a bosonic quantum field system 
at finite temperature, and at zero and non-zero chemical potential. For 
a simple spatial partition we derive the corresponding mutual informa- 
tion, a quantity that measures the total amount of information of one 
of the parts about the other. In order to find it, we first derive the 
geometric entropy corresponding to the specific partition and then we 
substract its extensive part which coincides with the thermal entropy of 
the system. In the case of non-zero chemical potential, we examine the 
influence of the underlying Bose-Einstein condensation on the behav- 
ior of the mutual information, and we find that its thermal derivative 
possesses a finite discontinuity at exactly the critical temperature. 
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1 Introduction 



The notion of the geometric entropy, in the framework of a quantum system, 
is quite old. One of the first calculations [1] was performed in the eighties for 
the case of a scalar field propagating in a black hole background. Some years 
later, a similar problem, in the framework of a quantum field theory, was ad- 
dressed by several authors [21 [3l HI [5] . Geometric entropy, generally speaking, 
is a measure of the information loss after cutting out a spatial region of the 
system. It caught attention because of its characteristic behavior: for a sys- 
tem in its ground state, it grows like the boundary surface of the excluded 
subregion, a property it shares with the black hole entropy. In fact, in the 
context of quantum field theory, pioneering work on the geometric entropy 
was driven in part by the suggested connection to the Bekenstein -Hawk- 
ing black hole entropy 0. From the very beginning, geometric entropy has 
been tightly related with the presence of spatial entanglement in a quantum 
system. Entanglement is a fundamental ingredient of quantum mechanics 
leading to strong correlations between subsystems. From the early days of 
quantum mechanics up until now, it has been playing an increasingly impor- 
tant role in understanding and controlling quantum systems. The interest in 
it has been renewed [7] after the developments of the quantum information 
science in which it is viewed as a resource in quantum information processing. 
Geometric entropy has been considered as a measure of spatial entanglement 
when the system under consideration is in a pure quantum state with a den- 
sity matrix of the form p = |^)(^'|. By defining an "in" and an "out" spatial 
region and tracing out the "in" degrees of freedom one obtains the reduced 
density matrix for the "out" region: pout = trin/5. The geometric entropy is 
then defined as the von Neumann entropy: Sout = — ti"yOout Inpout- 

When the system is in a thermal state, e.g. in a mixed state, the geo- 
metric entropy can be defined, following the von Neumann definition, in an 
analogous way. However, in this case, it does not have the same properties 
as the entanglement entropy in a pure state system, and it is no longer a 
good estimator of entanglement since it mixes correlations of difi'erent types 
[31[5l[8], from genuine quantum to thermal correlations. Since it measures 
the thermal information loss, geometric entropy becomes an extensive quan- 
tity at the limit of an infinite system, and loses the "area law" behavior that 
characterizes a pure state system. As an alternative probe for the amount 
of correlations between different parts of a system in the case of thermal 
states, the notion of the so-called "mutual information" [9] has been proposed, 
which, roughly speaking, eliminates the contribution of the extensive part of 
the thermal entropy from the geometric entropy and can be considered as 
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an upper limit for the entanglement entropy. 

In any case, geometric entropy has been considered as a convenient con- 
struction, playing the role of an order parameter, for the investigation of 
finite temperature conformal quantum field systems [5] and in the context 
of the AdS/CFT correspondence, aiming at the physics of strongly coupled 
Quark-Gluon Plasma, the weakly coupled deconfined phase of Yang-Mills 
theories or the phase structure of large QCD at a finite density |10| . 

In the current work, we examine the geometric entropy in a free bosonic 
quantum field theory at finite temperature and at zero and non-zero chemical 
potential. Having found it, we subtract its extensive part, that is, the part 
related to the amount of information that is lost due to the mixed nature 
of the system. The result, defined as the mutual information, quantifies the 
spatial correlations between the different parts of the system and exhibits 
the known area law behavior. 

The underlying reason for the present study is connected to the Bose- 
Einstein condensation that characterizes the system, which has been in the 
center of theoretical and experimental investigations during the last fifteen 
years after the production of the condensate in the laboratory Bose- 
Einstein condensation has the characteristics of a phase transition albeit, 
theoretically at least, it can take place in an ideal system |12 1 113 | [M j I15|. It 
is then natural to search for the interconnection between this phase transition 
and the spatial correlations in the Bose system. Our findings indicate that, 
indeed, the Bose-Einstein condensation influences the behavior of the mutual 
information: We find that its derivative with respect to the temperature, 
dIm/dT, has a finite discontinuity at the critical temperature both at the 
non-relativistic and the relativistic limit. Thus, we show how this phase 
transition leaves its fingerprint on a quantum informational quantity like 
mutual information. 

The paper is structured as follows: In section II we present the calculation 
of the geometric entropy at finite temperature and zero chemical potential 
and we define the resulting mutual information, setting the stage for our 
result goal which is the topic of the section that follows. This is section III, 
in which we apply our results in an environment with a finite charge density, 
we calculate the explicit form of the mutual information, and we derive the 
discontinuity of its temperature derivative. Finally, in the last section we 
summarize our findings. In Appendix A we present some technical details of 
our calculations. 
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2 Geometric entropy at finite temperature 



Our starting point is the thermal density matrix p = e /Z of a quantum 
field system and its Fock space representation: 

p[^',^] = ^{^'\e-^^\^), (1) 

where ^ denotes a single scalar field or a collection of fields. The matrix 
element ([T]) can be written as a functional integral: 

1 r - fdr fcPxCl^] 

p[$',$] = — - / V^T,x)e . (2) 



#(0,x)=$(x) 
$(/3,x)=$'(x) 

It worth noting that at the zero temperature limit, /3 — ?• oo, ([2]) is just 
an interpretation of the ground state density matrix [1]. It is then natural 
to expect that our result will reproduce, at this limit, the known [H IH [5] 
entanglement entropy. To derive the geometric entropy we follow the usual 
line of reasoning and we divide the D dimensional space on which our system 
is defined, into two regions A : {xi > 0,x±) and B : (xi < 0,x±). Tracing 
out the "in" region, and gluing along the axis xi > 0, n copies of the resulting 
reduced density matrix, we find [3l|4l[5]: 

In Zn the fields are defined on a D + 1 dimensional space M„ = Rd-i x C'n- 
The subspace Rd-i is an Euclidean space with metric ds^ = dx"^ + • • • + dx"]^ 
while Cn is a two dimensional Riemann space consisting of n sheets glued 
together along the positive xi axis. This n folded structure turns eventually 
[1] the (t, xi) plane into a flat cone with an angle deficit 5 = 27r(l — n) at 
the origin. Having found tr{pji)'^ the geometric entropy is defined through 
the relation: 



lim ^^^^^^ — - = -tv{pR In PR) = Sg = - (-^ - l] In Z„ 
n-fi n — 1 \on J 



■(4) 

n=l 
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For a free bosonic theory, the partition function can be deduced by fohowing 
standard steps: 

oo 

0+ 
oo 

X /^'-'■"'''0.e-<-«) (5) 

0+ 

where = 5^ + (?| + m?. Due to the locaUty of the action, the partition 
function in ([5]) is not expected to depend expUcitly on the details of the 
Riemann surface. Thus, in order to calculate the non-trivial trace appearing 
in (jSj, we start with the finite temperature propagator of a free particle in 
cartesian coordinates: 

A^^{x',x) = {a^\e-n-9l^\x)fs^. (6) 

The subscript /3,„ = /3n indicates the periodic boundary conditions imposed 
on the thermal Green's function. As it is obvious, they are dictated by the n 
folded structure of the Riemann space C„. The next step is to transfer the 
result onto a two dimensional cone with angle deficit27r(l — n): 

ds^ = dp^ + p'^n^dO^, < < 27r. (7) 

One can easily find that the free thermal propagator in ([6]) assumes the form 

m-. 



oo 

^ 4T 



e 

4ttT 

u=—oo 
oo 



{x'^-xiY + (x'f^-xo-vl3„Y 



1 2 
rrT ^ ' ^ ' 



U= — QO 



The above expression can be written in the conical metric ([7]) by making 
the replacements xq = psm{n6), xi = pcos{n6), and using the expansion 



e 

m=— oo 
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where J\-m\ are Bessel functions of the first kind. Thus, the thermal propa- 

n 

gator on the surface d?]) reads: 



1 



E 



e 4it e X 



i/,m,mi,m2 



p1 



(10) 



In the last expression the rotation T ^ it has been adopted in order to 
secure convergence of all our intermediate steps. 
Tracing out (|10p we find: 



A 



/3n 



-n-di) 



- E 

2it ^ 



._M 

e 4it i n 



P -(11) 



At this point we stress the fact that for n ^ 1, the trace over the conical 
metric ([7]), that is the integration over p, must be performed before the 
summations over m or mi. The relevant calculations can be facilitated by 
using the fact that we are only interested in the derivative of (jlip with respect 
to n: 



(8»-4ft.)„.i = E" "*" ^» 



E 



_(ifni 

e iit 



dppj 



.pL ._m . ( p^\ 

" J|m| I — j 





oo 



+ 



n=l 



71=1 



In obtaining the last expression we have used the identities: 



(12) 



(13) 



In Appendix A we prove that: 



oo 



pL ._m t f 

2itl n J,~ — 



2tJ it-)-T 



-4^ + l^U--)+0(^)(14) 
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and 



oo 



dppJl 



2t 



it^T 4-kT \V2 



(15) 



In the above equations we introduced an upper cutoff R in the p-integrals and 
we have written as V2 = vri?^ the volume of the two dimensional subspace. 
Substituting the first term in the rhs of (|14p into (jl2p and feeding with 
the result (|l'2p we find (see Appendix) that it leads to the logarithm of the 
partition function: 



vr — J T— V 



(4 



(16) 



Following the same steps for the first term in the rhs of (jlSp we can prove 
that it is connected to the thermal entropy of the system: 



1 VD-1V2 



dr, 



(47r) 2 



■0+ 



T 2 



_(pp_ 

(2^ 



ln(l 



(17) 



where oj"^ = + m? . The contribution to (fT2|) of the second term in the rhs 
of (1141) assumes the form: 



1 



1 



12 



(47r: 



dT 
T— 



-Tm- 



E 



e 4T 



Vd- 



d^p 



(27r)^wtanh(^)' 



(18) 



Collecting everything together and using ([5]) we get the geometric entropy: 

d^p 1 



o 



-Vd 



(27r)^u;tanh(^f 
d^p 



(2vr) 



D 



ln(l 



+ 



/3a; 



(19) 
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At the limit p — )• oo , tanh ^ VP — j _^ ^ consequently the fist in- 
tegral in (|19p diverges. The same divergence appears in the case of zero 
temperature: 

5,(/3 = oo) = |v„_,/^-^=L= 

After this observation we are led to write: 

= S,W = OO) + ^Vo-l I ^^^^l^ ^ + ^gthe.n.al. (21) 

Some comments are in order at this point. The first term in the last ex- 
pression represents the well known [H 131 [5] entanglement entropy at zero 
temperature. This is a quantity that diverges in the absence of an ultravio- 
let cutoff, while it grows like the boundary surface of the excluded subregion. 
This fact clearly indicates the existence of very strong quantum correlations 
between fields defined at neighboring points, a direct consequence of a lo- 
cal quantum field theory. A quantitative explanation of such a behavior 
can be traced back to the uncertainty relations. Even at zero temperature, 
the notion of a sharp, well defined, boundary surface is more classical than 
quantum. The divergences appearing in Sg{T = 0) are connected to the 
fact that in (I20p we integrate down to zero distance, driving to infinity the 
density of the reduced density matrix eigenvalues. The second term is finite 
and well-defined for m? > 0. It is also proportional to the boundary surface 
and it is an increasing function of the temperature. We can consider it as a 
measure of the number of degrees of freedom that have been excited on the 
boundary surface due to the non-zero temperature and, consequently, as a 
measure of the thermal correlations between the partitions. The last term 
is the thermal entropy of the subsystem, an obviously extensive quantity. 
Subtracting this term from the geometric entropy we are led to define the 
mutual information: 

IM = S,(fi = OO) + / (22) 

In general the mutual information is a measure of all correlations, thermal 
and quantum. We use the following definition: 

I{A:B) = S{pa) + S{pb)-S{pab)- (23) 
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For the case in hand the entropy of the combined system AB is just the 
total thermal entropy, S{pab) = •S'thcrmal- The entropies of each one of 
the two subsystems are equal due to the way we have divided our system. 
Moreover, each one of them contains a part which is one half of the total 
thermal entropy of the system. Thus, their extensive thermal contribution 
to the mutual information is equal to S{pab) and when substracting the 
latter, all contributions due to the thermal entropy will be eliminated. So, 
what (|22p represents is the mutual information of the system divided by 2, 
ImiP) ^ HA : B)/2. 



3 Mutual information and Bose-Einstein condensa- 
tion 

Almost all of the technical details needed for the current section have al- 
ready been exposed in the previous one. The basic difference of the analysis 
that follows, lies on the fact that we are now interested in charged scalar 
(non- interacting) fields. The field theoretical description will be based on 
complex fields while the introduction of a chemical potential (as a Lagrange 
multiplier) will ensure the conservation of the charge. In this framework, the 
partition function of the system assumes the form: 

Z{P)= J V(l)V(j)*expl - j dr J d^xC[(l),(j)*]\. (24) 

/3— periodic 

The Lagrangian entering the last expression can be written [131114) as follows: 
r ] = (a, - pf - dl + in"] 4>. (25) 

Following the same steps as in the previous section, we find: 



oo 



ln^„(/3) = ^j^Vb^. [ ^e-^-^trc^e-^^-^l+^^^o-/^^). (26) 

Once again we start from the free thermal propagator in Cartesian coordi- 
nates: 

1 °° f 1 

AisA^',x) = ^ ^ exp<^ __(f'_f)2 + 

1/=— oo I. 

+ [^ + f^)ixo-xo) — piyf3n> (27) 
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to arrive at the traced quantity that is relevant for the final calculation in 
(1261): 



trc„e-'^(-^^+2/^9o-M') 



ntrciC' 



-T(-9|+2/iao-M^) 



+ 



+ 



- - 

12 I n 



n 



(28) 



where the arrow underlines the fact that we have kept only the terms that 
are relevant for determining the mutual information. 
In Appendix A we show that: 



1 



(2^)^^ 1 tanh 



tanh 



2 



(29) 



Isolating the (diverging) zero temperature contribution we find: 



5,(/3 = oo) + ^Vd^i I 



d^p 



{2tt)^uj I e(^-/^)/3 - 1 



+ 



+ - 



1 



(30) 



For the system in hand the zero temperature entanglement entropy reads: 



Sgifi = oo) 



1 

6 (47r 



d^p 



D-1 
2 



D — 1 m^-x 
A^l 



(31) 



To reveal the physical content of our results we shall focus on the well- 
studied D = 2> case which hosts the Bose-Einstein condensation. As it is 
well-known [121 lT3l IT^ the quantitative realization of the phenomenon is 
different at the two opposite limits, the non-relativistic p <C and the 
ultra-relativistic one p ^ m^, as these are defined by the total charge density 
of the system. 

Beginning from the non-relativistic case, in which the charge density is 
very low and the anti-particle contribution can be omitted |14| . we rewrite 
(1301) in the form: 



tNR 



SgiP = oo) + 



TT V2 f d^p 



1 



6 m J (27r)3 



(32) 
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where we noted as i^nrH^) = fJ, — m < the non-relativistic chemical poten- 
tial. The integral appearing in (f32|) is the total density of particles occupying 
excited states: 



1 



\ 3/2 



m 



2vr/3 / 



E 



--NR 



n=l 



3/2^ 



(33) 



where Z]sir = e'^^^^ < 1. For temperatures below a certain critical value 
Tc we have n{Tc) = mn and the above quantity is a constant: 



, < T < Tc. (34) 



At exactly the critical temperature the number p4p becomes the con- 
served total particle density of the system: 



P 



27rm \ 



\ 3/2 



(35) 



As an immediate consequence we get for the mutual information: 

3/2 



tNR 



(36) 



Above the critical temperature the system passes to the gas phase in 
which all of the particles occupy excited states. The mutual information 
reads now: 



NR 



SgiT = 0) + ^^p, T>Tc. 
b m 



(37) 



Thus, the Bose-Einstein condensation and the relevant phase transition are 
reflected in a discontinuity of the derivative of the mutual information: 



(9T 



t=t; 



9T 



T=T7 



(38) 



When p we are approaching the ultra-relativistic limit, the critical 

temperature rises at relativistic high values Tc = (3|/9|/m)^/^ ^ m and the 



11 



behavior of the system changes. Below the critical temperature, one easily 
finds that [T2l[T3]: 



Sg{T = 0) + ^V2 I 



d^p ( I 1 

(27r)3w i e('^-/^)/^ - 1 ^ - 1 



(39) 



The integral that appears in ()38p and (j39p is not the charge density of the 
system and, consequently, is not a conserved quantity even for temperatures 
above the critical one. However, it is not hard to confirm [12| that at high 
temperatures T > Tc it behaves as following: 



d^p ( 1 
(27r3)a; I e('^-/^)/^ - 1 



+ 



1 



T 2 

(m 

6 2tt^ 



+ -i^(m2 
47r2 



g(ij+/i)/3 
2 



(40) 



where C = e"^^" /47r and 7^; is the Euler-Macheroni constant. 

As in the non-relativistic case, the derivative of the mutual information 
with respect to the temperature possesses a discontinuity that reflects the 
underlying phase transition: 



dirr 



dT 



dl„ 



t=t; 



dT 



t=t; 




(41) 



The last result completes our study for the influence of the Bose-Einstein 
condensation on the entropy of entanglement in an ideal Bose system at 
finite temperature and non-zero chemical potential. 



4 Conclusion 

In the current work we have performed two types of calculations and we 
have arrived at results with a clear physical content. First, we calculated 
the geometric entropy in an ideal Bose system at finite temperature and we 
confirmed the expected result: It combines the genuine quantum correlations 
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with the thermal fluctuations, and it becomes an extensive quantity for an 
infinite system. Due to the simplicity of the system under consideration, 
we were able to explicitly subtract from the geometric entropy its exten- 
sive component which coincides with the corresponding thermal entropy. In 
this way we defined the so-called mutual information, which grows like the 
surface that bounds the space region in which a system lives. The second 
calculation we performed refers to a Bose system at finite temperature and 
non-zero chemical potential. We found that, at the critical temperature, the 
temperature derivative of the mutual information exhibits a finite disconti- 
nuity, and we explicitly calculated it. This result connects the condensation 
that appears in an ideal quantum Bose system with the spatial correlations 
between two regions of the system. This connection was shown by using a 
purely informational tool namely, the quantum mutual information. 
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Appendix 

In this Appendix we shall prove those of the formulas appearing in the text 
for which summations over m or v must be performed. To begin with, let 
us discuss ()14p . The relevant integral diverges and calls for the introduction 
of a cutoff: 



To handle the last integral we make an intermediate step by introducing the 
following expression: 



oo 






(43) 
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which is also a regularized version (for a — t- 1"^) of the integral entering ([T 
Taking the limit a = 1 + e, e— t-O^ one easily finds that: 

We immediately see that the diverged part of the integral appears for n = 1. 
In this case the integration in (|4'2|) is trivial and we are led to the conclusion: 

1 vrii^ V2 , , 

Combining this identification with the finite part appearing in ()44p we get 
the confirmation of (jl4p . 

Our next concern is (|16p . Using the identities: 

(27r)^+i 

E e^^o^^ = y E 5(P0--.), -. = ^ (46) 

!^= — 00 ''^ fc= — 00 

we recast the integral appearing in Eq. (|16p into the form: 



e 4T = (47rr) 2 / -— ^e 



00 



' -e > e 4T 



2 



D + l I 13+3 ^ 

V 



(47r)— J T— 



^ 0+ 

/3 00 

To obtain the last result we wrote V2 = J dr J dxi — )■ f3L, we rescaled 

-00 

T — )• and we used the abbreviation = + . 

Performing the integral over T and neglecting an irrelevant (infinite) 
constant we get: 



00 



^g-T[(/3^)2 + (27rfc)2 



T 



-ln((/36j)2 + (27rA;)2). (48) 
0+ 

The summation over k is standard 1131: 



^^ln{if3ojf + i27Tkf) = l/3oj + lnil-e-f^'^). (49) 



2 

k 
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The last result proves ()16p . 

Following the same line of reasoning we can prove (jlSp . Using, once 
again, the identities (j46p we rewrite the relevant integral in the form: 



1 Vd- 



12 



D-l 

(47r) 2 



TT 



dT 

D + l 

T— 



-Tm' 



E 



e 4T 



+ a; 



2 • 



The summation is easily performed: 

1 13 



(50) 



(51) 



Combining (j50p and (I5ip we immediately obtain (llSp of the text. 

Our next concern is Eq. ([T5|). The relevant integral diverges and the 
introduction of a cutoff is necessary. To this end let us discuss the integral: 

oo 



E 



-I- 2lo tanh(/3a;) 



^ / dppe ^Jln\ 



2t 



-P 



8*2 



(i//3i?)^ 



(52) 



which can be considered (at the limit R — t- oo) as a regularized version of the 
integral appearing in (fTSl) . Note that the divergence in (I52p is independent 
of n and, contrary to (j44p . there is no finite part for n 7^ 1. This completes 
the proof of (115p . 

To prove (17) it is enough to follow the road we followed to arrive at 
(j49p . Beginning from the relation: 



1 VD-1V2 f dT „y„2 ^ 

O , ^ D+l / D+3 ^ / J 

2(4^)— Vr— V 



e 4T 



4^ y (27r)^ 2^ 



^p-T[(/3„a;)2+{27rfc)2] 



(53) 



0^ 



we only have to perform a differentiation with respect to n to arrive at the 
result indicated in (fTT 



d^p r 1 



n 

_dPp_ 
{2tt 



(2vr) 



D 



In 1 



-Pu)n — In ^1 



-[Sum 



n=l 



(54) 
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The last relation we have to prove is (p9|) of the text. We begin by using 
the Poisson summation formula to find that: 



oo 



k 



^g-TK+iM)^ (55) 



/5 . 

With the help of this result we get the mutual information: 
2tt f d^p ^1 



3p 'J {2tt)D ^ {oJk + ill? + 00^' 

The sum in the last expression can be easily performed if we rewrite it in 
the form: 

y \ 

^ (wfc + inY + up- 

l a;-^ ^p 1 1 a; + ^p 1 

2 uj ^2 + _ ^,)2 + 2 uj Z^u;l + {u + f,)^- ^ > 

Using for each term the formula (|5ip we get the result indicated in (I29p . 

It would be useful to compare our result indicated in (jl9p with the cor- 
responding result derived in the framework of a two dimensional conformal 
scalar field theory |5j with central charge c = 1/2. This can be done by 
identifying the ultraviolet cutoff A with the inverse lattice spacing 1/a and 
the mass m with the inverse finite size of excluded interval l/l (that is, the 
infrared cutoff). Given that our result is valid at the limit L?"m?' — t- oo, the 
comparison is meaningful only for /3/a — t- oo or ///3 — t- oo. Applying (jl9p for 
D = 1 we get the result: 

f i In f a") ^ 1 In f l") if /3 ^ oo 

6 m/3 ^ 6/3 



if /3 
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